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Comment on "Quasiparticle Spectra around a 
Single Vortex in a d-wave Superconductor" 

In a recent Letter jjj Morita, Kohmoto and Maki 
(MKM) analyzed the structure of quasiparticle states 
around a single vortex in a d x -i_ y 2 superconductor using 
an approximate version of Bogoliubov - de Gennes (BdG) 
theory with a model d-wave gap function. Their prin- 
cipal result is the existence of a bound state within the 
core region at finite energy with full rotational symmetry, 
which they assert explains the recent scanning tunneling 
microscopy results on YBa2Cu307_a single crystals ||. 
Here we argue that the approximate BdG equations used 
by MKM are fundamentally inadequate for the descrip- 
tion of a d-wave vortex and that the obtained circular 
symmetry of the local density of states (LDOS) is an 
unphysical artifact of this approximation. 

BdG equations for a d-wave superconductor differ from 
the conventional ones by the non-locality of the off- 
diagonal term f dx'A(x, x')u(x'). MKM approximate 
this term by a local term A(k;R)u(R), and claim that 
this is accurate to 0(l/fcp£), quoting the derivation given 
by Bruder jjj. A key step in Bruder's calculation is fac- 
toring out the rapid kp oscillation in w(x) by introducing 
a slowly varying envelope function Vk(x): 


u(x) 


Wk(x). 


(1) 


v(x) is then substituted in the off-diagonal matrix el- 
ement and after transforming to the center of mass 
frame both A and v are Taylor-expanded in the rela- 
tive coordinate. The leading term in the expansion is 
that used by MKM and corrections are indeed formally 
0(1/ kp£)- Since one cannot expand the rapidly varying 
function u(x) directly, this manipulation rests heavily on 
the ansatz (1). An implicit assumption made in (1) is 
that the wave-vector k is a good quantum number. This 
is approximately correct in the simple planar scattering 
situations considered by Bruder 0. However, k is cer- 
tainly not a good quantum number near the vortex core; 
the appropriate ansatz in this situation would involve a 
sum over k in Eq. (1). One then obtains the leading 
off-diagonal term of the form £ £ e lfcF£R A(k; R)u k (R), 
indicating that all the plane wave states are mixed to- 
gether by the non-trivial structure of the gap. This in 
turn results in coupling of the various angular momen- 
tum channels and the loss of rotational symmetry. MKM 
improperly neglect the summation over k in the ansatz 
(1) for the wave functions. 

Inadequacy of MKM's approximation becomes even 
more transparent when analyzing its implications for the 
physical quantities, such as the LDOS. A brief examina- 
tion of the (standard) BdG equations for a d-wave vor- 
tex reveals that both the quasiparticle amplitudes and 
the self-consistent gap function will necessarily exhibit a 
four-fold anisotropy. This qualitative conclusion is cor- 
roborated by the numerical results within the Eilcnberger 
H and the lattice BdG || formalisms. Anisotropic 
wave- functions are also found near a strongly scattering 


non-magnetic impurity ||. These anisotropics are the 
salient feature of unconventional superconductivity re- 
sulting from an intricate interplay between the center of 
mass and the relative coordinate in the off-diagonal part 
of the BdG Hamiltonian. The fact that such anisotropics 
have not yet been observed experimentally J2| is indeed 
puzzling and apparently inconsistent with the above the- 
oretical work. However, reconciliation between the the- 
ory and the experiment cannot be brought about by sim- 
ply ignoring the symmetry breaking interactions as is 
done by MKM. 

A brief inspection of MKM Eqs. (1) and (2) reveals 
that they in fact describe a set of decoupled s-wave vor- 
tices with the gap functions of various amplitudes (but 
identical profiles) scaled by the factor cos[20(kYj. Each of 
these vortices possesses a set of bound states Q| with the 
energy spacing A 2 (k; oo)/Ep, contributing sharp peaks 
in the LDOS. The finding of MKM of the broad peak in 
the LDOS in the vicinity of the core is a trivial conse- 
quence of superimposing the peaks corresponding to the 
lowest bound states belonging to these s-wave vortices. 
The full rotational symmetry of the LDOS in MKM's re- 
sult also directly follows from this observation. Clearly, 
such calculation does not resolve the interesting question 
of the existence and the structure of the quasiparticle 
bound states in d-wave vortices. We expect that the na- 
ture of such states will be qualitatively different from the 
conventional picture of Caroli, de Gennes and Matricon 
0. However, this structure can emerge only from a solu- 
tion, numerical or analytical, which preserves complexity 
and richness of the original BdG theory. 
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